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Abstract 

We show that comparison of the detailed calculations of the synchrotron radiation of expanding magnetized 
clouds (plasmoids) with the temporal and spectral evolution of radio flares observed from the recently discov¬ 
ered galactic “microquasar” GRS 1915+105 imposes strong constraints on the model parameters of relativistic 
plasmoids, in particular, on the absolute values and time-dependence of the magnetic field, speed of expansion 
of radio clouds, the rates of injection of relativistic electrons into and their energy-dependent escape from the 
clouds, etc. The data of radio monitoring of the pair of ejecta not only enable unambiguous determination 
of the aspect angle and the speeds of propagation of both components, but also contain an important infor¬ 
mation about the energy source for acceleration of the electrons, in particular, may distinguish between the 
scenarios of bow-shock powered and magnetized relativistic wind powered plasmoids. Within the framework 
of the synchrotron-self-Compton model we also calculate the fluxes of synchrotron hard X-rays and GeV/TeV 
inverse Compton 7-rays, and discuss constraints on the parameters of GRS 1915+105 which could be provided 
by positive detection or flux upper limits of these energetic photons. 

1 Introduction 

The galactic black hole (BH) candidate source GRS 1915+105, discovered in 1992 as a hard X-ray 
transient (Castro-Tirado et al. 1994), attracts an extreme interest since the detection of superluminal 
radiojets in this source by Mirabel & Rodriguez (1994; hereafter MR94). Together with the second 
galactic superluminal jet source GRO J1655-40 ( Tingay et al. 1995; Hjellming & Rupen, 1995), they 
constitute a new class of objects called microquasars (MR94), representing a scaled down analogs of 
QSOs (or generally AGNs) in two principal aspects: both microquasars and quasars are most probably 
powered by accretion onto BHs, although of essentially different mass scales, and both populations 
produce relativistic jets. Being, however, much closer to us than AGNs, the microquasars enable radio 
monitoring of not only approaching but also receding ejecta (importantly, in short timescales), which 
makes them unique cosmic laboratories for study of the phenomenon of relativistic jets (MR94). 

Here we summarize the basic results of our recent study (Atoyan & Aharonian, 1997; AA97) of 
the prominent March/April 1994 radio flare of GRS 1915+105, when the approaching (south) and 
receding (north) ejecta moving in opposite directions with a speed /? ~ 0.92 at an aspect angle 6 — 70° 
have been resolved, and superluminal speed 1.25 c (for d = 12.5kpc) has been discovered (MR94). 

2 General scenario 

Comprehensive discussion of radio flares observed from GRS 1915+105, which most probably are 
connected with ejection of radio emitting material (i.e. clouds, containing relativistic electrons and 
magnetic fields frozen in the co-moving plasma), is given in MR94, Rodriguez et al. (1995), Foster et 
al. (1996). Here we would emphasize two important features of time evolution of the flares: 

(a) a rapid, less than a day, rise time of the fluxes, S u oc zz~ Qr , revealing transition from optically thick 
(a r < 0) to thin (a T > 0) emission, which is a typical signature of an expanding radio source; 

(b) the spectral index a r ~ 0.5 at the stages of flare maximum later on steepens to a r ~ 1. 


1 


The study of these features results in conclusive information on the intrinsic parameters of the 
radio source. Indeed, modeling the ejecta (in its rest frame) as a spherical source of a radius R and 
magnetic field B , with a power-law distribution of relativistic electrons IV( 7 ) oc 7 -2 extending to 
7 > 10 4 , we can express R and B in terms of S'* = S 10 GHz/500 mJy and cloud’s opacity r v = Rx k u 
(k u is the synchrotron absorption coefficient) at a frequency v as (AA97): 

B ~ 0.15 S^rj-^S^ 17 ^ 17 ^ 7 G , (1) 

R = 4.6 x 10 14 r 28 / 1 V 1/17 ^ A V/^gh / 17 cm , ( 2 ) 

where vghz = vf 1 GHz, 5 = \J\ — /3 2 /(I — (3cos 9) ~ 0.57 is the Doppler factor, and 7 = w e /wB is 
the ratio of the electron to magnetic field energy densities. For the outburst of 19 March 1994, the 
fluxes measured at instant to — 4.8 days after ejection (the time of VLA observations on March 24.6), 
corresponded to S* ~ 1. Therefore, as it follows from Eq.(2), for optical transparency 77 < 1 apparent 
down to v ~ 1.4 GHz, one needs R> 8 x 10 14 cm. To reach that size during t' 0 = 5 to ~ 2.7 days, the 
radio clouds should expand in their rest frame with a speed u exp > 0.1 c. This agrees with the speed 
of expansion deduced directly from observations, and suggests that the lack of detection of red-shifted 
optical lines might be due to large Doppler-broadening of the lines (Mirabel et al. 1997). 

For T 1 . 4 GHz < 1 Eq.(l) results in B < O. 27 ” 4 / 1 ' G at times t — to- For the equipartition (7 = 1) 
magnetic field B eq ~ 0.2 G the synchrotron cooling time of electrons emitting at frequencies 12 ~ 
10 GHz exceeds few years. Therefore the synchrotron losses cannot be responsible for the observed 
steepening of radio spectra to a r — 1, which requires a steepening of the electron spectrum IV( 7 , i) 
to an index a e — 3. We suggest a model (AA97) which attributes the steepening of IV( 7 , t) to (i) 
continuous injection of relativistic electrons into the radio cloud, with the rate < 3 ( 7 +) oc 7 ~ 2 q(t), and 
(ii) energy-dependent escape T esc (j, t) oc 7 -1 of the electrons from the cloud. 


3 The March/April radio flare of GRS 1915+105 


For calculations of the nonthermal radiation expected during radio flares, we have to find temporal 
evolution of the energy distribution of electrons IV(y, t) in an expanding magnetized cloud. Generally, 
the equation for IV( 7 , t), including the therm for stochastic in-situ acceleration, represents the partial 
differential equation of the second order (e.g. see Kardashev 1962). However, in the cases when the 
region of particle acceleration can be distinguished from the region of bulk emission, the acceleration 
terms can be effectively replaced by injection of accelerated particles into the radiation production 
region, which results in a well known differential equation of the first order: 
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For the parameters characterizing the radio clouds in GRS 1915+105, the energy losses P = 
— (dy/dt) of electrons in an expanding magnetized medium are mainly due to adiabatic losses, and 
at very high energies also the radiative losses. Assuming the temporal evolution of magnetic field 
in the form of B(t) = Bo(R/Ro)~ m , where R = R(t) and Ro = R(t-o) , we have P = u exp 7 /R + 
P 2"/ 2 /R 2m (p ‘2 = const). Assuming u exp = const, we have found (AA97) analytical solution to Eq.(3). 
This solution can be used in numerical calculations for an arbitrary time dependence of v exp (t), if 
approximating the latter by the mean v\ in small time intervals A+ 

Allowing for a (natural) deceleration of cloud’s expansion at later stages, we suppose v exp = 
vq /(1 + t'/t exp ) k , where t! is in the rest frame of the ejecta. The injection rate is taken as <3(7+0 oc 
exp(— 7 / 7 c ) q(t') , which implies that the shape of energy spectrum (power law with exponential 
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Fig.l a (left): Top panel - The speeds of approaching (solid line) and receding (dashed line) ejecta calculated 
from Eq.(4) for different ratios of S a /S r . In the symmetric case, /3 a = /3 r (= 0.920), the aspect angle is 9 = 69.6°, 
resulting in S a /S r = 12.8; Bottom panel -Evolution of the radio fluxes from the south (S) and north (N) radio 
components expected in the case of identical ejecta with /3 a = /3 r (dot-dashed lines), and slightly asymmetrical 
ones, with S a /So = 7 (solid lines) and 9 (dashed lines). Note that the total fluxes coincide. 

Fig.l b (right): The time evolution of the radio fluxes and spectral indexes at 8.42 GHz, expected in the case of 
synchronous increase (by factor of 2 ) of the injection rate of relativistic electrons into the south and north radio 
clouds during At < lday around equal intrinsic times t' = 9days (solid curves). The dashed curves show the 
spectral evolution without an ‘afterimpulse’. The dotted curves correspond to the case when the ‘afterimpulse’ 
is in the form of an increase of magnetic field by 40%, while the injection of electrons proceeds smoothly. 
To distinguish the flux curves, somewhat different shapes of the ‘magnetic’ and ‘electronic’ afterimpulses are 
supposed. The principal difference between these 2 cases is in different time evolution of a r . 


cutoff) does not change in time. The term q(t') = (1 + £ , /ii n j) -p defines the time profile of the electron 
injection rate. For the escape time we suppose energy dependence in the form r esc (7, t') oc R~ u j~ A 
and take into account that r esc > R/c. 

Model parameters 

Comparison of the numerical calculations of synchrotron radiation of radio clouds with the fluxes 
detected during March/April flare (MR94; Rodriguez et al. 1995) prove that initial speed of expansion 
vo ~ (0.1 — 0.2) c is needed, depending on the magnetic field Bo at instant to = 4.8 days after ejection. 
For the equipartition magnetic field, B eq ~ 0.2 G, a typical power of injection of relativistic electrons 
L e ~ (0.5 — 1) x 10 38 erg/s is required. However, L e should be much higher if Bo < B eq , e.g., 
L e > 10 39 erg/s for Bo — 0.05 G. 

The observed rate of decline of the flare requires a decrease of B(R ) oc R~ m with m — 1. The 
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power law index a ; n j ~ 2 follows from the observation of a r — 0.5 on March 24, when the escape losses 
have not yet significantly modified the electron spectrum N( r y,t). The observed steepening of the 
radio spectrum at ~ 10 GHz from a r ~ 0.5 to a r ~ 0.84 on April 16 requires, in our model, an escape 
of the electrons with A ~ 1, and allows for some variation of the parameter u ~ 1.5. Importantly, 
the steepening of radio spectra can be explained only if the cloud’s expansion after f exp ~ (1 — 3) days 
gradually decelerates, with k ~ (0.7 — 1), and if about the same time (namely, tj n j ~ (1 — 2) t ex p ) the 
injection of new particles declines with p ~ 2k. Otherwise the escape of radio electrons from rapidly 
expanding cloud would not be fast enough to steepen N( 7 , t) as needed. Remarkably, an injection 
rate of this kind can be easily provided, in particular, if one would assume that it is proportional to 
the solid angle of the cloud as observed from the central source: q(t') —► %(t / ) = {R(t') / v^t'] 2 . We 
prefer to call qb(t') as the ‘beam injection’ case, and further on use this type of injection in order to 
reduce the number of free model-parameters. 

Brightness ratio as a measure of South/North asymmetry 

To explain the observed (MR94) brightness ratio of the approaching to receding ejecta S a /S r = 8+1 
(instead of S a /S T > 12 expected for identical plasmoids), we suppose that the twin ejecta are similar 
but not identical (for another possibility see Bodo & Ghisellini, 1996), allowing for some asymmetry 
between the pair of ejecta in the speeds of propagation (/3 a 7 ^ /3 r ), aspect angles (0 T / 180° — 0 a ) 
or intrinsic luminosities ( S' a 7 ^ S' T ). In particular, the degree of South/North asymmetry needed is 
smallest if mainly the speeds /3 a and f3 r would be different. In this case for the ratio S a /S r at equal 
intrinsic times we have 

Sa = = / 1 + A- cos 9 \ 3+ ° r 

5b UJ vrj V1 Pa COS 6) ' 1 j 

Due to strong dependence of Eq.(4) on the ratio of Lorentz-factors T r /T a , a deviation of the ratio 
S a /S r from the one expected for identical ejecta can be explained by a small difference between /3 a and 
A (see top panel in Fig.la). For example, solution of Eq.(4) together with the equations describing the 
apparent angular separation of the ejecta, results in A. = 0.928, A = 0.905, 6 = 70.3° for S a /S r = 7. 
In this case the fluxes of both components seem to agree with the ones measured by MR94 (see Fig. la, 
bottom panel). Note that for the given parameters of GRS 1915+105, the ratio S a /S r = 7 at equal 
intrinsic times corresponds to 5 a / S r ~ 8 measured at equal angular separations. 

Synchronous afterimpulses far away from the central source ? 

The reported accuracy of VLA flux measurements in MR94 is ~ 5 %. Meanwhile some of the 
data points in Fig. la, in particular the fluxes on 21-st day (April 9) for the south, and 28-th day 
(April 16) for the north components, exceed the calculated fluxes by 5 %, therefore they still need a 
better explanation. We would like to believe that this discrepancy could be connected with significant 
increase of the fluxes between April 4 and 5 (t = 16 — 17days) detected by the Nancay telescope. 
It is seen from Table 1 and Fig.l in Rodriguez et al (1995), that during 24 hours between these 2 
days the fluxes at both 1.41 GHz and 3.3 GHz have suddenly increased by ~ 30%. Although VLA 
was not observing GRS 1915+105 at that days, the ‘echo’ of that event could be present in the flux 
detected by VLA from the approaching component on April 9. Also, one should expect significant 
delay between the times of VLA observations of that event from the south and north radio clouds, 
if it was due to powerful two-sided afterimpulse from the central source which has reached the two 
clouds at equal intrinsic times f a = t' r . In Fig. lb we show the fluxes calculated for practically the 
same model parameters as in Fig.la (for 5 a /5 r = 7), but assuming that there was an additional 
short (At 1 < 1 day) impulse of injection of relativistic electrons into both clouds around intrinsic times 
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t' = 9 days after ejection event. Taking into account different Doppler factors of counter ejecta, this 
intrinsic time corresponds to apparent times t a = 16.6 days and t r = 27.5 days for the approaching and 
receding components, respectively. It is seen from Fig. lb that agreement with the measured fluxes 
now becomes better. Note that the last two data points (on April 30) cannot be explained by another 
afterimpulse, since they coincide in time. The excess fluxes on that day are most probably connected 
with the second ejection event occured, as noted in MR94, around April 23. 

If such an interpretation of the data is not just an artefact, but corresponds to reality, implications 
for the physics of jets may be very important: it could mean that both clouds continue to be energized 
by the central source being even far away from it (presumably, through a relativistic wind in the 
jet region, and perhaps also the wind termination reverse shock on the back side of the clouds), or 
otherwise we have to rely upon a mere coincidence of equal intrinsic times (as well as amplitudes) of 
additional injection of relativistic particles from the bow shocks ahead of two counter ejecta. 

4 Predictions for nonthermal high energy radiation 

GRS 1915+105 is a very strong transient X-ray source, with peak luminosities exceeding 10 39 erg/s. 
This radiation most probably originates in the thermal plasma near the central compact object. In 
addition, nonthermal X-rays of synchrotron origin could be produced in the relativistic ejecta, if the 
spectrum of accelerated electrons there would extend beyond TeV energies as it is the case in blazars 
(e.g. see Urry &; Padovani 1995). Although from the point of view of required energetics in jets it is 
obvious that the X-radiation of the jet generally cannot be responsible for the bulk of observed X-ray 
fluxes, the synchrotron component, however, may become visible at high energies, say Ex > 100 keV, 
where the thermal emission is typically suppressed. 

In Fig.2 we present the time evolution of the synchrotron and IC radiation fluxes in the range from 
radio to very high energy 7 -rays expected in the case of cloud’s magnetic field Bq = 0.05 G, calculated 
for the injection spectrum of electrons with exponential cutoff around E c = 30 TeV. It is seen that 
during the first several hours after ejection (or somewhat more, depending on the actual time needed 
for formation/acceleration of the expanding ejecta), when the cloud is still opaque for synchrotron 
emission at GHz frequencies, the synchrotron fluxes may dominate over the extrapolation of thermal 
component at energies > 100 keV. This may result in significant flattening of overal spectrum at 
Ex ~ 100 keV. Interestingly, hard tails of X-ray spectra perhaps are a common feature of the galactic 
BH source population, among which is GRS 1915+105 (Harmon et al. 1994). Although an existence 
of such feature in the spectrum of GRS 1915+105 could be seen only marginally (see e.g. Sazonov et 
al. 1994), in case of the second microquasar GRO J1655-40 the X-ray spectra clearly extend up to 
several hundred keV (see Harmon et al. 1995). 

The second signature of acceleration of TeV electrons in GRS 1915+105 could be detection of IC 7 - 
rays. The calculations in the framework of synchrotron-self-Compton models, based on the parameters 
defined from radio observations, show that during the strong flares we may expect detectable fluxes 
of 7 -rays , first of all at TeV energies, provided that the magnetic held in the radio clouds does not 
exceed 0.1 G. In Fig.2 we show the fluxes of IC 7 -rays calculated for parameters of the March/April 
1994 hare at 3 different times after the outburst. We can see that the 7 -ray fluxes above several 
hundred GeV during first few hours of the hare are on the level of 3 x 10 ” 11 erg/cm 2 s, which is about 
of the VHE 7 -ray hux of the Crab Nebula. After few days the hux drops to the level of 0.1 Crab which 
is still detectable by current Imaging Cherenkov telescopes in Northern hemisphere (CAT, HEGRA, 
Whipple). Afterwards the source becomes invisible. 

Assumption of lower magnetic held results in higher (approximately oc R -2 ) 7 -ray fluxes. However, 
as far as a magnetic held significantly less than 0.05 G would imply injection luminosity of electrons 
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Fig.2 : Broad band synchrotron (heavy lines) and IC (thin lines) radiation fluxes expected during a flare from 
GRS 1915+105 at t = 0.1 day (solid lines) t = lday (dashed lines) and t = 10 days (dot dashed lines) after 
relativistic ejection event. Data points correspond to the level of maximum radio fluxes observed at 8.4 GHz 
(the cross) and 234 GHz (diamond) (from Rodriguez et al. 1995). The star is for the flux of IR jet reported by 
Sams et al. (1996), and the hatched region shows the level of hard X-ray fluxes typically detected during the 
X-ray flares of GRS 1915+105 in the energy range > 20keV (Harmon et al. 1997). Magnetic field Bq = 0.05 G 
is supposed, which implies L e = 3 x 10 39 erg/s. 


L e > 10 10 erg/s, it is difficult to expect y-ray fluxes essentially exceeding the fluxes shown in Fig.2. 
On the other hand, in the case of magnetic field on the level of equipartition B ~ 0.2 G or higher, 
the IC fluxes dramatically decrease. Therefore either positive detection or upper limits of VHE y-ray 
fluxes, being combined with hard X-ray observations above 100 keV, could provide robust constraints 
on the magnetic field in the ejecta (or the luminosity L e ) and/or efficiency of acceleration of electrons 
beyond TeV energies. 

Information about IC y-rays can be obtained also in EGRET energy domain, i.e. at E > 100 MeV. 
However, it should be noted that even during the first few hours of the flare, the fluxes in this energy 
range do not exceed 10~ 6 ph/cm 2 s. This implies that with the effective detection area of the EGRET 
~ 10 3 cm 2 , only several photons could be detected during the first day of the outburst. Unfortunately, 
since the IC fluxes significantly drop with expansion of the cloud, observations of the flare after that 
time could not noticeably improve the photon statistics. 
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